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Abstract. We realize on an Atom-Chip a practical, experimentally undemanding, 
tomographic reconstruction algorithm relying on the time-resolved measurements of 
the atomic population distribution among atomic internal states. More specifically, we 
estimate both the state density matrix as well as the dephasing noise present in our 
system by assuming complete knowledge of the Hamiltonian evolution. The proposed 
scheme is based on routinely performed measurements and established experimental 
procedures, hence providing a simplified methodology for quantum technological 
applications. 
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The estimation of quantum states by using measured data is crucial to verify the 
quality of any quantum device. To fully determine a quantum state, i.e. to perform 
a quantum state tomography, one needs to accumulate enough data to compute 
the expectation values of an informationally complete set of observables [1]. The 
availability of a complete set of measurements to be implemented by the experimenter 
is not straightforward and in general quantum state reconstruction is carried out by 
complicated set-ups that have to be robust against noise and decoherence sources, in 
order not to limit the accuracy of the reconstruction Eli- 

The standard technique used in quantum optics for the full reconstruction of 
quantum states is coherent homodyne detection |3] , which has recently been extended to 
ultracold atoms in [5]. Homodyne detection requires the use of a, not always available, 
local oscillator held that acts as a phase reference for the state under reconstruction 
[ 0 ]. For atomic internal states simpler interferometric techniques can be used to map 
the relative phases of internal components onto the level populations [7j. Cold atomic 
systems and degenerate quantum gases are indeed unique tools for quantum simulations 
[6] and precision measurements of atom characteristics beyond the classical limit [8]. 
However, their applications outside laboratory depend critically on the simplihcation 
and downsizing of bulky cold-atom setups. Indeed, an interesting technological 
development is given by the possibility of integrating cold atoms with nanostructures [9] . 
An important invention in this direction is the realization of a Bose-Einstein Condensate 
(BEC) in microscopic magnetic traps based on the micro-electronics technology, yielding 
the so-called Atom-Chip m HI]- For example, experiments based on Atom-Chips 
have demonstrated nonlinear interferometers with the sensitivity beyond the standard 
quantum limit [8], nonclassical interferometry with motional states [12], and quantum 
Zeno dynamics na. 

This work presents a practical, experimentally undemanding, tomography protocol, 
that relies only on the time-resolved measurements of the atomic population distribution 
among atomic internal states. This protocol allows the reconstruction of a, not 
necessarily pure, state of a n-level quantum system, where the coherence elements are 
unknown and usually challenging to be measured. The idea is simple: the state to be 
reconstructed evolves in time and the population distribution is measured at different 
times; the same dynamics is numerically simulated starting from a randomly chosen 
initial state; we run an optimization protocol that minimizes the difference between 
the simulated and measured data in such a way that the optimal solution provides the 
tomographic reconstruction of the initial state. Our procedure is based on the complete 
knowledge of the system evolution. Notice that, by applying the proposed scheme to a 
complete set of known states, one is in principle able to reconstruct the system dynamics, 
leading to a procedure similar to the quantum process tomography na. 
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The experimental apparatus is based on a microscopic magnetic trap, or atom chip 
m [16], where we bring an atomic sample of ®^Rb to quantum degeneracy by forced 
evaporative cooling. Most of the structures and wires necessary to the magnetic 
trapping, the forced evaporation and successive manipulation of the atoms are all 
embedded in the atom chip, making this device versatile and experimentally easy to 
use HT]. Our BEG has typically 8 10^ atoms in the low held seeking hyperhne state 
|F = 2,mF = 2), at a critical temperature of 0.5 /iK and is 300 /rm away from the 
chip surface. The experiment described in this work is performed 0.7 ms after the BEG 
release from trap to guarantee a homogeneous magnetic bias held and strongly reduce 
the ehects of atomic collisions. In this way the most relevant source of noise on the 
evolution turns out to be the instability of the environmental magnetic held. 

We consider a u-level system represented by the hve fold F = 2 hyperhne ground 
level of ®^Rb. In the presence of a magnetic bias held the degeneracy between the 
magnetic sublevels is lifted according to the Breit-Rabi formula. Using two external 
Helmholtz coils, we arbitrarily set the magnetic held to 6.179 G [TH]. To drive the atomic 
dynamics we apply a radio-frequency held (RE) oscillating at urf = 27r 4.323 MHz [T9] 
using a wire structure embedded in the chip. Thanks to the relative proximity between 
the BEG and the emitting wire with ~ 10 mW of RE power we excite magnetic dipole 
transitions between the atomic levels at Rabi frequencies up to 200 kHz. 

To record the number of atoms in each of the mp states of the F = 2 hyperhne 
state we use a Stern-Gerlach method. After the state manipulation has been performed, 
in addition to the homogeneous bias held, we apply an inhomogeneous magnetic held 
along the quantization axis for 10 ms. This causes the diherent mp states to spatially 
separate. After an expansion time of 23 ms a standard absorption imaging sequence is 
performed. The recorded atomic population in each mp state is normalized to the total 
number of observed atoms. 

It is important to stress that each measurement sequence completely destroys the 
system. This means that each sampling is associated to a diherent experimental cycle 
of production and manipulation of the BEG and we have to rest on the assumption that 
both the state preparation and the tomographic procedure always yield the same result. 
This assumption is a posteriori conhrmed by the agreement between the prepared and 
reconstructed state. We sample every 5 oscillation of the RE held, corresponding to 
having a data point every ~ 1.16 /is. We repeat 5 times each measurement to obtain 
the mean values of the relative atomic populations = Pi{tj) and standard deviations 
(Tjj = (Tiitj) for each sublevel i. An example is represented in Figj^ 

3. State reconstruction 

The post-processing analysis is formulated in the following way: suppose that we want to 
estimate the initial quantum (not necessarily pure) state of a n-level system, described by 
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Figure 1. Time evolution of an arbitrary atomic state. An initial unknown state 
evolves under the effect of the Hamiltonian H and the five atomic populations (points) 
are recorded every 1.16 /rs, i.e. every 5 oscillation of the RF field applied to the atomic 
sample. The superimposed line is the theoretical evolution with the estimated initial 
condition po minimizing the error e(po) between experimental and theoretical data. 


the density operator po = pit = 0). To do that, we use the a priori known Hamiltonian 
operator H describing its dynamics and the measurements of the subsystem populations 
at different evolution times Piit). In other words, we measure the observables aja*, hence 
obtaining the expectation values Tr[p(t)ajaj] = Pi{t), where Oj and a\ are, respectively, 
the annihilation and creation operators for each sublevel i, and p{t) is the time-evolution 
of the unknown state po that we want to estimate. In the case of fully coherent 
Hamiltonian evolution, one has p{t) = U{t)poU{t)^ with U(t) = exp[—iHt] being the 
unitary evolution operator. If the system is, instead, subjected to a noisy evolution, as 
it is the case in most experimental situation, one has p{t) = ^t{po) with <I)t being the 
so-called quantum map or quantum channel or quantum operation mapping po into p{t) 

m- 

In our case the Hamiltonian, written in the rotating wave approximation, is 
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where the state basis is chosen to go from mp = +2 to mp = —2, the RF held Rabi 
frequency is H = 27r 60 kHz, the detunings Si are dehned by = 3 kHz and ^2 = 11 kHz. 
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Figure 2. Pictorial representation of the theoretically expected density matrix pin 
(left) and the reconstructed one po (right), with the real (upper) and the imaginary 
(lower) components. The Uhlmann fidelity J-{po,Pin) between po and p™ is = 0.98. 


In order to include the unavoidable presence of dephasing noise, mainly originated 
in our experiment by the presence of magnetic held huctuations superimposed on the 
bias held, we add in our model a Lindblad super-operator term £, acting on the density 
matrix p as C{p) = P} + 2 ajaj p ajoj], which randomizes the phase 

of each sublevel j with a homogenous rate 7 . We hnally obtain the density matrix 
evolution as 

j^p[t:) = -'pH,p(t)]+C(p(t)). (2) 

To reconstruct the state one has — 1 unknown independent real parameters 
dehning the nx n density matrix, with the constraint of leading to a physical state, i.e. 
a semi-dehnite positive, hermitian matrix p of unitary trace. This problem can be also 
formulated in the language of semidehnite programming izu. It corresponds to hnding 
the initial state po minimizing the diherence (e.g., mean squared error) between the 
measured populations pj(fj) and the quantities Tr[<l>t(po)alhi] for each observed time step 
tj, hence 11 Tr[<hq (po)alai] - Pi{tj)\\ with ^t{po) mapping po into p{t) being 

the solution of the master equation in ([^, and || ■ || being any mathematical norm. 
Here, we have implemented the minimization algorithm using the Subplex variant of 
the Nelder-Mead method [ 22 ] • In particular, to take into account also the experimental 
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Figure 3. Behaviour of 1 — Pin) versus the tomography time (used to minimize 
the error e(po)), where J-{po, Pin) is the Uhlmann hdelity between po and po is the 
state reconstructed from the experimental data, and pi„ the expected one. 


uncertainties, we minimize, with respect to po) Ih® following weighted mean squared 
error function 



( 3 ) 


where Pij(po) = (Po)dlai] and cuij = 


4. Results 

To test our scheme, we have prepared, by applying known Hamiltonian evolutions, a 
set of states to be reconstructed. For example, in Fig. we reconstruct a state that 
is obtained by applying a 7r/2 pulse [23] to the \F = 2,mF = 2) state. We report 
the experimentally recorded population evolutions in a 16 ps-long time window. We 
then compare these results with the theoretical evolutions of the diagonal elements (i.e. 
populations) of the reconstructed state po- For this reconstruction, based on 16 averaged 
observations of the hve-level populations, the computed error is e(po) ~ 2 x 10“®, 
corresponding to an Uhlmann fidelity [2l] of F{po, Pin) = 0.98 between the reconstructed 
density matrix po and the prepared one pi„, that are pictorially represented in Figj^ 
Let us point out that no a priori knowledge of the initial state has been used for the 
tomographic reconstruction. However, this information has been exploited to calculate 
the Uhlmann hdelity. 

Moreover, we have faithfully reconstructed, with errors below 3 x 10“®, other states 
prepared by using the Hamiltonian in Eq. ([^ and randomly varying over time the 
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detunings and 62 - These low errors correspond to fidelities higher than 0.95. 

To check the reconstruction error convergence with respect to number of collected 
data we apply the optimization algorithm to the population distributions in different 
time windows T, and computed the quantity 1 — -T(po,Pm) in each case. The results 
in Figj^show that, when T is comparable with the system natural evolution timescale 
(T ~ 27r/ff), the reconstruction accuracy is already satisfactory, with the error quickly 
saturating to its minimum value. 

Furthermore, this technique can be also exploited to get further information on the 
system evolution, e.g. estimate the amount of dephasing noise in the system dynamics 
resulting from its coupling to the external environment. 



Figure 4. Error e{po) as a function of the tomography time Ti (reconstruction window 
[ 0 , Ti\) and the dephasing rate 7 , whose optimal values ^opt minimizing the error e are 
shown in the dashed line. 

More specifically, we have reconstructed the initial density matrix measuring the 
population evolution within different time windows [0, Tj], with Tj G [10, 100] /is, 
but neglecting the presence of dephasing, i.e. 7 = 0 . Then, for any reconstruction, 
we have found the optimal value of dephasing rate 7 opi G [0, 750] Hz minimizing the 
error e between the theoretical and experimental data - see Fig. It turns out that 
the optimal dephasing rate increases with the tomography time Tj, hence selecting the 
noise spectrum components larger than 1/Tj, as expected. In other words, another 
important application of our procedure is to extract information on the presence of 
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external noise and then indirectly on the strength of the coupling between the system 
and the environment. We expect this to allow the characterization of the type of external 
noise, with the possibility of identifying the presence of temporal correlations. 

5. Conclusions 

To summarize, we have experimentally demonstrated a tomographic reconstruction 
algorithm that relies on data collected during the evolution of an unknown quantum 
state, assuming a complete knowledge of the system Hamiltonian. The advantages 
of this protocol are the simplicity of the post-processing procedure and the use of a 
quite conventional absorption imaging technique. Furthermore, we have shown the 
convergence of the protocol even using a small amount of collected data, compared 
to standard tomographic technique. Finally, we have also estimated the rate of the 
dephasing noise present in our system dynamics by repeating this procedure for longer 
tomography time windows and minimizing the reconstruction error. 

Moving another step further, a promising application may be represented by 
the possibility of characterizing the noise itself, for instance its spatial and temporal 
correlations, by investigating the behaviour of the state reconstruction error in terms 
of different noise models. The proposed scheme therefore realizes quantum state 
tomography but could readily be modified to perform quantum process tomography 
by assuming complete knowledge of the input states, hence providing a very feasible 
and useful tool for several quantum technological applications. 
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